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Abstract 



A partially first-order form of the characteristic formulation is introduced 
to control the accuracy in the computation of gravitational waveforms pro- 
duced by highly distorted single black hole spacetimes. Our approach is to 
reduce the system of equations to first-order differential form on the angular 
derivatives, while retaining the proven radial and time integration schemes of 
the standard characteristic formulation. This results in significantly improved 
accuracy over the standard mixed-order approach in the extremely nonlinear 
post-merger regime of binary black hole collisions. 
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I. INTRODUCTION 



One of the most impressive achievements of the characteristic approach is its demon- 
strated ability to carry out long-term stable evolutions |] of generic single black hole space- 
times. Work is under way to extend these numerical simulations to the post-merger regime 
of binary black hole coalescence, by first computing the gravitational radiation emitted dur- 
ing a white hole fission 0, and then extending these results to compute the gravitational 
waveforms emitted during the post-merger phase of a binary black hole collision P|. 

The present work is an outgrowth of one of these projects, having been motivated by 
numerical difficulties that we have encountered with the present implementation of the 
characteristic code JiJ. In the course of our simulations of a white-hole fission ||, we have 
noticed an angular oscillation mode in some of the metric variables of the characteristic 
formulation, similar in appearance to a red-black || decoupling. Numerical experiments 
reveal that although the area that it affects can be somewhat reduced, its amplitude does 
not converge to zero with increasing grid resolution for practical grid sizes in the rapidly 
changing environment near the time of merger of a binary black hole. We had expected to 
see convergence in this regime, given the measured second order convergence of the code on 
less rapidly changing solutions, even at fairly low resolutions @J|. 

This effect is clearly of numerical nature, and would not have been observed in earlier 
work, as it is only triggered when, in addition to extreme nonlinearities in the metric func- 
tions (such as those found in the post-merger phase of the head-on collision of two black 
holes, or in the time- reversed picture, during the fission of a single black-hole), a fairly large 
degree of asymmetry and of time dependence in the metric functions is present. Moreover, 
it appears that even when some of these conditions are satisfied, the dissipation built into 
the algorithm used for the evolution the conformal metric of the spheres |7j can, in certain 
circumstances, be sufficient to suppress the spurious oscillations, provided the boundary 
data is specified analytically |§. 

In the systems that we have previously considered, such as Robinson- Trautman space- 
times, a perturbed Kerr black hole || and the collapse of matter into a black hole 0, our 
choice of the gauge conditions at the boundary effectively precluded fulfilling the second 
condition above. In all cases, the boundary data was given analytically, so no extraneous 
effects could be observed. 

This oscillatory mode is not an unstable mode of the finite difference approximation in 
the usual sense of exponential growth of the error. This is confirmed by tests where random 
initial and/or boundary data || is prescribed, and the code run for many crossing times. 
But it is worth pointing out that while those tests are useful to spot exponentially growing 
modes, they would not detect a non-exponentially growing oscillatory mode such as the one 
we address here. 

We observe, in fact, an inaccuracy in the computation of the metric variable W, which 
does not converge to zero, and gets triggered in situations where there are large angular 
variations and rapid time changes in the metric variables, specially in /3, Eq. (|T|). Because 
(3 is related to the time slicing of the boundary (and in the black hole case, this boundary 
is the event horizon 7i + ), (3 will indeed change substantially near the time of coalescence 
of a binary black hole horizon (or, in the time-reversed description actually used by our 
numerical code, near the time of fission of a white hole). 
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The PITT null code runs stably in the usual sense (without exponentially growing modes) 
in the configurations that we have explored. The effect of the angular noise is nevertheless 
important for our present purposes, since the angular error it introduces in the metric 
functions is most pronounced near null infinity, and this in turn spoils the accuracy of the 
waveforms computed. 

We have found that the difficulty lies with the characteristic evolution equations them- 
selves, and it is not a boundary-induced effect. Although the boundary (horizon) data is 
smooth and, by construction, satisfies the consistency conditions fllO| , the evolution and 



hypersurface equations themselves introduce the above mentioned oscillation mode, which 
is shown clearly in Fig. [1|. We have found clear evidence that the origin of the difficulty lies 
with those terms in the equations that contain second angular derivatives. These type of 
terms are present both in the evolution equation for the conformal metric of the surfaces 
at constant luminosity distance r, and in the hypersurface equation for the metric function 
W itself. Most notably, when this terms are suppressed, so is the angular oscillation mode 
in W. The most troublesome terms are those containing second angular derivatives of the 
metric variable /3, of the form SB/3 and S 2 /?. 

Second angular derivatives of the conformal metric Hab, Eq. (0), also appear in the 
hypersurface equation for W , in the form of 8 2 and 56 operators acting on quantities defined 
in terms of the conformal metric. However, these terms do not seem to lead to numerical 
difficulties as pronounced as those caused by the terms involving second angular derivatives 
of /3, although the reason behind this is not entirely clear at the moment. It is possible that 
in the present regime the contribution from these terms is not as important. 

Our rationale for implementing the characteristic code in mixed first and second order 
form has been that doing so is entirely straightforward, and leads to an accurate and stable 
discretization with the least number of variables. The argument of using the least variables 
is compelling; however one must consider whether relaxing this requirement might not lead 
to better overall numerical behavior. We will show that this is indeed the case. 

Rigorous stability arguments can be given for the linearized system in the axisym- 
metric case, with the analysis extending unmodified to the full three-dimensional problem as 
shown in || . Later work |7| has shown that special care must be taken in the nonlinear case, 
as some unstable modes are present that are not revealed by the Von-Neuman analysis of 
the linear problem. One must keep in mind however that a stable discretization of a model 
problem as in [7| may not address other issues such as those raised here. The numerical 
analysis for mixed order systems is not yet at the level of sophistication achieved for first 
differential order systems. Long-term stable evolutions |] have been achieved, and this is a 
consequence of the remarkably good choice of variables and coordinates inherent to Bondi's 
system of equations, as well as of the care taken in implementing them numerically. 

We have attempted to solve the problem at hand by the expeditious method of intro- 
ducing further dissipation in the equations, in the spirit of . Regrettably, these numerical 
remedies do not seem to have the desired effect, apparently because they do not go to the 
root of the problem. 

We consider here a different solution, which is analytic in nature, but motivated by a 
numerical application, specifically the need to attain controlled accuracy in the computation 
of gravitational waveforms out of highly distorted single black hole systems. Our approach is 
to remove the offending second angular derivatives from our system of equations altogether, 
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by introducing additional variables, in terms of which the enlarged system can be written 
in first differential order form in the angular coordinates. We thus implement a partial 
reduction to a first differential order system, leaving the radial part of the evolution equations 
in second order form, and allowing for mixed radial and angular derivatives, both in the 
original equations and in the auxiliary equations introduced. One goal is to preserve, where 
appropriate, the radial and time integration schemes developed for the PITT null code, 
whose long-term stability (in the usual sense of not having exponentially growing modes) 
has been amply demonstrated 0. At the same time, we perform the necessary modification 
to the system of equations to make the PITT null code useful for the problem at hand. 

The plan of the present paper is as follows: in Sec. [TI| we briefly review the characteristic 
formulation in its standard second-order form. This material has appeared in Refs. HH, 
and we include it here in a condensed form, as needed for the derivation of the partially 
reduced system. Our new results are contained in Sec. |T| and onward, where we present a 
partial reduction of the characteristic formulation, with all second order angular derivatives 
eliminated with the aid of a (minimal) set of additional variables. In Sec. [IV], we review 
the discretization strategy for the system of equations introduced in Sec. [TTI[ Finally, we 
present an application of the partially reduced system in Sec. [V] to the computation of the 
spacetime exterior to a white hole horizon, comparing our results with those obtained with 
the standard null cone formulation. 



II. NULL CONE FORMULATION IN SECOND ORDER FORM 



The material in this section has appeared in Refs. P|| which should be consulted for a 
detailed derivation. Here we present the main equations (including all the nonlinear terms), 
in a slightly more compact form than that of Ref. H, for ease of reference and to provide a 



starting point for the approach developed in Sec. |n. 

We use coordinates based upon a family of outgoing null hyper surf aces, letting u label 
these hyp ersur faces, x A (A = 2,3), label the null rays, and r be a surface area coordinate, 
such that in the x a = (u, r, x A ) coordinates, the metric takes the Bondi-Sachs form ||T~2|.[T3 



ds 2 = - 



r J 



,W I 1 + ) _ r 2 h AB U A U B 



du 2 - 2e 2p dudr 



2?" h^sU dudx + r h A sdx dx 



(1) 



S A 



where W is related to the Bondi-Sachs variable V by V = r + W; and where h AB \iBc 
and det{h,AB) = det(q A B), with q AB a unit sphere metric, given in terms of a complex dyad q A 
satisfying q A q A = 0, q A q A = 2, q A = q AB q B , with q AB q B c = $c and Iab = \ (qAqB + UIb)- 
We represent tensors on the sphere by spin- weighted variables Jl4| . The conformal metric 
h A s, which satisfies the condition det(h A B) = det(gAs), is represented by the complex 



1 + JJ. The 
i|,fl. Angular 



function J = \h A Bq A q h \ and by the real function K = ^h A sq A q B , where K 2 
metric functions U A are similarly encoded in the complex function U = q A U A 
derivatives are expressed in terms of S and 5 operators, for details, see Ref. 

The equations for the null cone formulation, in second order form, are obtained directly 
from the appropriate components of the Ricci tensor MM . These equations form a hierarchy, 



U 
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splitting into hypersurface equations, which involve only derivatives on the null cone, i.e. 
R rr = provides an equation for /3 r in terms of J, while R r Al A = gives U rr in terms of J 
and /3, and the trace R,ABh AB = yields W >r in terms of J, (3 and U, respectively. 
The main equations MM consist of the hypersurface equations 



P,r = g (JrJr ~ K) > 

(r 2 Q) >r = -r 2 (8j + 3^), r + 2r 4 3 (r~ 2 /3) r 

+ r 2 (jl - fr)(3K r + 9J r ) + 3(JJ r ) + 9(Jif >r ) - J r 8K 

£/, = — (#Q - JQ) , 

^e 2/3 ft - 1 - e^clcV 3 + V 2 (r 4 (3C7 + 8C/)) r 



+ e 



2,; 



1 



(1 - K){mi3 + 3/33/3) + -( J(3/3) 2 + J(5^) 2 



i (p/3{dK - 3 J) + d(3(QK - 8 J)) + i ( j8 2 /3 + J8 2 /3 



— e 



-2/3 ' 



-{2KU r U r + JU 2 r + JU 2 ). 



and the evolution equation for J, obtained from RabQ A( 1 B , 
2(rJ) iW .-(r- 1 V(rJ), r ) r = 

_ r -i ( r 2 5f/ ) ^ + 2r - 1 e 2 / 3 (3 2 /3 + (5/3) 2 ) - (r^W) r J + Nj. 

where 

K = 2K- mK + -(8 2 J + 3 2 J) + -777(8 J9 J - 8J3J). 

2 V ' AK ' 



(2) 



(3) 
(4) 



(5) 



(6) 



(7) 



and we have used the intermediate variable Q = r 2 e 2/3 (J£/ r + KU, r ) . Nj denotes the 
non-linear terms 

Nj = N n + N J2 + N J3 + N M + N J5 + N J6 + N J7 + -{P 1 + P 2 + P 3 + Pa) (8) 

r 



where 



e 213 / - - - - - 

iVn = — (K(5J8/3 + - 3J8/3) + J(3 J3/3 - 25^3/3) - J3J3/3 



N J2 = -- (3J(rC/ r + 2U) + BJ(rU r + 2Z7) 
N J3 = (1 - K){rW. r + 28*7) - J{rW. r + 23C7), 
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N M = r -e- 2 ^{K 2 Ul + 2JKU r U r + J 2 U 2 r ), 
N J5 = -^J r (W + W), 

N M = rQ(C73J + U3J)(JJ r - JJ. r ) 

+(JK, r - KJ ir )UdJ - U(QJ >r - 2KdKJ :r + 2JdKK, r ) 
-U(BJ r - KdJJ t r + JdJK tr ^J , 

N J7 = r(J tT K - JK, r ) (U(BJ - QK) + U(BK - 5 J) 

+k{w - w) + ( jw - jgc/)) , 

Pi = r 2 (^£(J r K - JK, r ) + ^(J r K - JK r )) - SV/3 t r, 

p 2 = e 2(3 (^- 2K(mp + 8/35/3) - (ppm + QftM) 

+ (j(8 2 /? + (Bp) 2 ) + J(3 2 /3 + (3/?) 2 )) + (8 J8/3 + 8 JS/3) J , 
^3 = ~ ((r9C/, r + 28C7) + (r9tT r + 29C/)) , 

4 

P 4 = - r —e~ w {2KU )T U tr + J£/ 2 + JE/J). (9) 

In Refs. [||,|| we have made a partial reduction to a first order system by introducing the 
auxiliary variable Q, which plays the role of U r . We sometimes refer to this system as first 
order (in time) because only the first (retarded) time derivative of J appears. See Refs. (||| 
for details of the discretization of the above system of equations. 



III. NULL CONE FORMULATION IN PARTIALLY REDUCED FORM 

A reduction to first order form of the system of equations in Ref. |J, in the quasi- 



spherical approximation, was considered by Frittelli and Lehner [JL5]. In their work, the main 
variables are Jab = h A B — Iab, Ua, ft and W = V — r + 2m, which measure deviation from 
Schwarzschild. The auxiliary variables they introduce to reduce the system to first order 
form are the following radial and angular derivatives of the metric functions, E>a = ft, a, 

M A bc = Jab,c, Q A = and P AB = Jab,t- 

Here we do not work with tensor components directly, but instead, following Refs. H|J, 
we represent tensors in terms of spin-weighted variables. For example, instead of writing the 
equations in terms of P>a, we work with B = q A ft,A = 8/3, and express its angular derivatives 
in terms of 8 and 8 operators. Second angular derivatives of ft can in turn be expressed in 
terms of <3B and 8P>. 

Eliminating all second angular derivatives in the hypersurface and evolution equations 
necessitates the introduction of exactly three complex variables. Since by symmetry there 
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are only three independent components of the conformal metric \iab-, there are then only six 
possible angular derivatives h,AB,Ci which are encoded in the three complex quantities 

^=dJ=h AB ,cq A q B q C , (10) 

v =dJ=~h A B,cq A q B f, (11) 

k = QK= h A B,cq A q B q C , (12) 

with spin weight (3), (1) and (1), respectively. The first order system for the quasi-spherical 
case considered in [ff3J introduces six new variables just for the angular derivatives of the 
conformal metric because no use is made there of the determinant condition. The con- 
formal metric h A B satisfies the condition det(/iAe) — det(gAs), it contains only two 
independent pieces of information, which are given in our notation by the complex function 
J = \liABq A q B ■ The remaining dyad component, given by the real function K = \h,ABq A q h '', 
is completely determined by the relation K 2 = 1 + J J, which is a consequence of the deter- 
minant condition. 

Therefore only two of the three complex quantities in Eqs. dlTD,(12) can be given in- 



dependently, with the remaining one fixed by the determinant condition. Our choice is to 
introduce the additional variables v = 8 J and k = (5K. By inspection of Eqs. (0)-©, we 
note that to put the system in first order form in the angular derivatives, we need only 
introduce one more auxiliary variable, B = 3/3. 

The new variables are initialized at the boundary with B = 3/3, v = 3J and k = dK 
respectively, and the consistency conditions 

B, r = 3/3 r> (13) 

",r = 3^,r, (14) 

fc r = m >r (is) 

propagate them to the interior, being the additional hypersurface equations for the new 
variables. 

The complete system of equations consists of Eqs. (7)-(12) of Ref. ||, together with 
Eq. flT5|), i.e., the evolution equation for J 



2 ( rJ ) ur ~ ( r ~ ly ( rJ ) r) = ~ K ( r ®Ur + 2 3f/) + -e w (m + B 2 ) 

V ' / ,r T 

- (rW, r + W^jJ + J H + JP U , (16) 
and the hierarchy of hypersurface equations 

V,r = cVr (17) 

k :r = QK, r (18) 

Ar = I (Jj,r - K) ( 19 ) 

B, r = 3Ar (20) 
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(r 2 Q), 



+ 



2K 2 



— K(k r + u r ) + vJ r + J9J r + z/K r + Jfc r — J r /c 

+ 2r 2 E r -4r5 



(j r - j 2 j r ) + gj(j r - J 2 J r ) 



Ef P = e 2/3 (KQ - JQ) 



l W) <r = K je 2/3 ipB + BB) + J (95 + B 2 ) + (v - k)Bj 

-l + 2rW + y 9C/ r - e~ w r —U, r {KU, r + JU, r ) I 



where 



ft = ft ( 2iT + 9 (v - k) + — (|9J| 



2 ) 



(21) 
(22) 

(23) 
(24) 



In Eqs. (pT6|) — (p4|) we have expressed the second angular derivatives in terms of 9 operators 
acting on B, v and k, replacing also 9 and 9 operators acting upon J, K and (3 with the 
corresponding auxiliary variables whenever possible. 

As in |J, we have regularized the W >T equation by setting W = r 2 W, and as in Appendix 
A of 0, we split the terms which vanish in the quasi-spherical approximation |J into two 
parts, one which contains only hypersurface derivatives, 

J H = — ( - KdJB + [Kv + (K 2 - 1)9 J - 2Kk) B 

(2k - v)B - 2K(QB + BB) + 2ft [{v — k)B + J (clB + 5 2 )] 



+ J 



+ r -e- 2 P( (KU r + JUrf 



U r (KU r + JU r ) 



- 2 [ p ( rU ,r + 2U ) + ® J ( rU ,r + 2C/)J + Ji3[9 (rU r + 2U)} - r J r ft[9C/] 
+ r (C/9J + C/i/) [Jj r ] - r (UQJr + t/^r) 

- 2r ( JK r - KJ, r ) (^[Uk] + i$s[KW - JW]) 

- 8 J ( 1 + rW) P tT , (25) 



and another, where we have explicitly isolated the only non-linear term where (retarded) 
time derivatives of J appear, 



2r 



Pu=Y^ [ J ,u(JrK-JK :T .)]. 



(26) 



The procedure just described eliminates all second angular derivatives from the hyper- 
surface and evolution equations. In computing the Bondi News, we must also evaluate the 
conformal factor u, which relates the metric on the x a = (u, r, x ) coordinates to the metric 
on a Bondi frame at X + ||]. (The relevant expressions are given in Appendix B of Ref. |4[] and 
we will not repeat them here.) Inspection of Eqns. (B1)-(B6) of || reveals that the second 
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angular derivatives of the conformal factor lo enter in the calculation as well. To remove 
these derivatives, we introduce the auxiliary variable W = <5u. Since u> is defined solely on 
Z + , we use the consistency relation W )U = Su; in to propagate W along the generators of J + , 
initializing it with Wq = <Suj , where u = u(x A ,u = u ). 



IV. NUMERICAL INTEGRATION 



We use the compactified coordinate x = r/(R + r), such that x = 1 at future null 
infinity, X + . Thus, the relevant radial derivatives can be written as d r — (1 — x) 2 /Rd x , 
rd r = x(l — x)d x and r 2 d r = Rx 2 d x , all of which are regular at x — 1. As in [[J, we carry 
out the radial integration with the right-hand side of Eqs. (|2l| ) evaluated at midpoints of 
the radial grid, 



Qi (Xi(l - Xi) + Ax) - <5,;_i {Xi-i(l - Xi-l) - Ax) 



(r 2 Q), r 



Ax, 



(27) 



and the same treatment is applied to Eq. (^). Note that the hypersurface and evolution 
equations are manifestly regular at 1 + when expressed in the variable x, except for the terms 
—ArB in Eq. (21) and 2rQU in Eq. (|23|) respectively, which have an extra factor of r. This 
factor, which would make them singular at X + , is effectively canceled by the corresponding 
factor of 1/r which appears in the right-hand-side of Eq. (^7|), hence the equations are regular 
for r > 0. Note also the limiting form of Eq. ([H]) at X + , Q = —2B, which is useful in the 
discretization of £7 r . For details on the numerical implementation of the eth operators, see 
Ref. 1 14]. The radial integration algorithm is explained in detail in Refs. ||TTl,|9|,|4lj7|1 and we 
will not cover it here. A departure from [[| is that we use a 3-step iterative Crank-Nicholson 
scheme |16| to ensure stability of the time evolution For the time integration of the 
conformal factor we use a combination of second-order Runge-Kutta and mid-point rule || 
integration schemes. 



V. WAVEFORMS FROM A FISSIONING WHITE HOLE 



We have implemented the system of equations presented in Sec. |TJ, and used it in 
the calculation of the waveforms emitted by a fissioning white hole. The details of that 
calculation are beyond the scope of this work and will appear elsewhere |2||| . The spacetime 
exterior to a white hole horizon is computed, by solving a double-null initial-boundary value 
problem. Boundary data is given on a white hole horizon, and initial data is specified on an 
outgoing null surface emanating from the horizon. 

A procedure for the construction of the complete boundary data, i.e. the intrinsic and 
extrinsic geometry of the white hole horizon, and a description of how to obtain the boundary 
data necessary to compute the exterior space-time via characteristic evolution can be found 
in Ref. [IIJ. Some regularization of the equations is needed to ensure we deal only with 



fields which are regular in the entire exterior spacetime ||. 

By reducing the system of equations presented in Refs. |I0|J2]] into the form presented in 
Sec. |HI| , with the addition of the auxiliary variables B = 5/?, v = SJ and k = <5K, we are 
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able to remove the angular oscillation clearly present in the metric function W when this is 
computed with the standard approach. 

We take as our benchmark case the most nonlinear (highest eccentricity) case considered 



in Ref. [fTOj] , with an eccentricity parameter of e = 10 . This corresponds to a case where 
the white hole horizon pinches off before the expansion of the outward null rays goes to 
zero, a necessary condition for a Bondi evolution forward in time throughout the pre-fission 
period. 

Figures [I] and |2] show the metric function W at X + on the north stereographic patch for 
both approaches. There is a clearly visible high frequency angular mode present in the edge 
of the patch in the calculation performed with the standard approach, as seen in Fig. [I]. 
This oscillatory mode is clearly absent from Fig. |2|. 

Using both approaches (standard and reduced), we have computed the waveforms, which 
correspond to the radiation emitted by a highly distorted white hole. In the time- reversed 
scenario, this would correspond to the radiation incident from X~, and the white-hole hori- 
zon, to the post-merger phase of a black hole horizon formed during a binary black hole 
collision. How to use this information to control the amount of incoming radiation and to 
gain insight into the radiation emitted during the post-merger phase of a black hole collision 



will be discussed elsewhere ]T8||2|,|5|] . 

Our motivation in introducing the partially reduced system was to control the numerical 
difficulties encountered in the computation of the radiation, as encoded in the Bondi news. 
The quality of the numerical waveforms is greatly improved in the partially reduced system. 
This is readily apparent comparing Figs. |3| and ^, which display the real part of the Bondi 
News, at a fixed time, for the standard and the reduced form of the equations. (In the plots, 
the region of the patch below the equator has been masked out.) While the evolution of 
the exterior spacetime is stable with the standard approach, the numerical noise swamps 
the signal. The noise is more pronounced near the edges of the stereographic patch, and 
it propagates to the interior. Figure ^ shows that, with the partially reduced system, the 
Bondi News at I + is perfectly smooth. 



CONCLUSIONS 

We have introduced a system of equations which takes the traditional null cone formu- 
lation and casts it as a system of equations which is of first differential order in the angular 
variables. 

The effectiveness of this approach is clearly demonstrated by comparing the metric func- 
tions at null infinity, and specially the waveforms at 1 + when computed via this approach 
versus those obtained with the traditional characteristic approach, in which second-order 
angular derivatives are present. 

We have thus shown that the partially reduced system of equations is highly effective 
in suppressing numerical errors, which otherwise would adversely affect the calculations of 
waveforms in highly nonlinear scenarios. 

We are currently using this approach to compute gravitational waveforms emitted during 
a white hole fission and during the post-merger phase of a binary black hole collision , 
and we are studying their application to the characteristic problem with matter sources. 
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FIGURES 




FIG. 1. The metric function W at null infinity, displaying the angular oscillation mode, a 
consequence of the second angular derivatives in the original system of equations. While the 
code is not unstable in the usual sense, the accuracy of the computation cannot be improved by 
increasing the grid resolution. 
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FIG. 2. The variable W at null infinity, when the second angular derivatives are removed by 
introducing the auxiliary variables B = S/3, v = 6 J and k = dK. The figure shows that the angular 
oscillation mode has been eliminated. 
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Re[N] 




FIG. 3. The real part of the Bondi news, 5ft [-/V], displays also a pronounced angular oscillation 
mode. The part of the north stereographic patch below the equator has been masked out. 
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FIG. 4. The real part of the Bondi news, R[iV] as computed from the metric functions obtained 
from the numerical evolution of the modified system of equations of Sec. III. Comparison with 
Fig. m makes it clear that the angular oscillation mode is no longer present. 
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